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We describe the occurence and properties of liquid crystal phases showing two dimensional splay and bend 
distortions which are stabilised by flexoelectric interactions. These phases are characterised by regions of locally 
double splayed order separated by topological defects and are thus highly analogous to the blue phases of 
cholesteric liquid crystals. We present a mean field analysis based upon the Landau-de Gennes Q-tensor theory 
and construct a phase diagram for flexoelectric structures using analytic and numerical results. We stress the 
similarities and discrepancies between the cholesteric and flexoelectric cases. 
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Elastic distortions in the nematic phase of liquid crystals 
can be categorised into three types: splay, twist and bend. Sta- 
ble phases with two dimensional twist distortions have been 
well characterised and observed experimentally. Here we pre- 
dict the occurence of structures with two dimensional splay- 
bend distortions. 

When nematics are doped with chiral molecules they can 
show stable phases with a natural twist known as cholesterics, 
shown in Fig. QI a )- In general cholesteric phases comprise a 
o«e-dimensional helix. This is because topological constraints 
mean that it is not possible to construct a state with helical 
ordering in two dimensions without introducing defects, or 
disclination lines, into the structure. However local regions of 
double twist are possible. Fig. |2ja) shows a simple, two di- 
mensional example of a director field where regions of double 
twist are separated by topological defects. If the free energy 
advantage of the double twist regions offsets the disadvantage 
of the disclinations phases of this type will be stable. 

Such states, the so-called blue phases of cholesteric liq- 
uid crystals, have indeed been observed [1]. The disclination 
structure is, however, more complicated than that of the sim- 
ple example in Fig. [2 a). The disclinations form textures with 
cubic symmetry and the stable phases have space groups 8 ~ 
and O 2 . The blue phases are stable at the isotropic-cholesteric 
phase boundary as here the magnitude of the order is small and 
hence the free energy penalty associated with disclinations is 
low. They are stable only over a narrow temperature range 
~1 K, although recently this has been extended to ~ 60 K by 
adding polymers or bimesogenic molecules Jj, □[ . 

One might therefore ask whether analogous behaviour 
could be observed in a liquid crystal which shows splay and 
bend, rather than twist, distortions. In 1969 Meyer showed 
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FIG. 1: Schematic representation of elastic distortions of nematic 
liquid crystals, showing the director field projected into the plane of 
the page, (a) Twist in the cholesteric phase of chiral nematics. (b) 
Splay and bend in a nematic with flexoelectric interactions. 
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FIG. 2: Schematic examples of two dimensional director field con- 
figurations, (a) Structure displaying regions of local double twist, 
characteristic of cholesteric blue phases, (b) Analogous structure 
possessing double splay distortions. 



that the one dimensional splay-bend distortion of the direc- 
tor field, shown in Fig. |TJb), could result from flexoelectric 
coupling to an external field (J] . Here we extend these results 
to show that near the isotropic-nematic transition two dimen- 
sional splay-bend structures (e.g., Fig. |2|b)) can be stable. 
We shall term such phases flexoelectric blue phases. Similar 
director field configurations, but stabilised by a saddle splay 
term in the free energy, have been reported in J5[] . 

The coupling of a liquid crystal to an external field occurs 
in two principal forms. The most usual is dielectric, where 
the field couples directly to the orientational order parame- 
ter. This coupling is quadratic in the field strength and the 
molecules tend to align parallel or perpendicular to the field, 
depending on the sign of the dielectric anisotropy. Flexoelec- 
tricity, which we shall deal with here, occurs because of the 
elastic properties of liquid crystals. If the liquid crystal is 
composed of pear-shaped molecules an elastic distortion can 
lead to a spontaneous polarisation. Conversely, if a polari- 
sation is induced by an external electric field then an elastic 
distortion results. This is termed flexoelectric coupling, and it 
gives a response linear in the field strength. 

Our conclusions are based upon Landau-de Gennes mean- 
field theory. We present analytic and numerical results ex- 
ploring the phase diagram of the flexoelectric blue phases, 
stressing the analogy, and differences, between these and the 
cholesteric blue phases. 

At a macroscopic level, the order parameter Q may be taken 
to be a traceless, symmetric, second rank tensor which is re- 



2 



lated to the anisotropic part of the dielectric tensor jfj. Using 
a tensor allows both the magnitude and the direction of the 
order to be recorded. The direction of the order, described 
by a vector n called the director, is defined as the eigenvector 
corresponding to the maximal eigenvalue of Q. We consider 
a bulk sample with periodic boundary conditions and assume 
that the liquid crystal has zero dielectric anisotropy and the 
electric field is uniform throughout the sample. 

The equilibrium thermodynamics of a nematic liquid crys- 
tal can be described by the Landau-de Gennes free energy [ 6], 
which we supplement by a familiar term for chiral coupling 
and an additional fiexoelectric coupling 
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Here, r is a reduced temperature, L% is an elastic constant, qo 
determines the pitch in the cholesteric phase, Ef is a fiexoelec- 
tric coupling constant, E is the electric field and V is the vol- 
ume of the domain, Q. There are four fiexoelectric couplings 
up to second order in Q, three of which may be written as a 
total divergence and converted to a surface term using Stokes' 
theorem. Since we neglect surface terms in this work only the 
term quoted need be retained. For simplicity, we adopt the one 
elastic constant approximation where the magnitude of twist, 
splay and bend are taken to be equal. 

Since chirality and fiexoelectricity induce structure in the 
fluid it is useful to analyse potential stable states by introduc- 
ing a Fourier decomposition of the Q-tensor 



Q(r): 
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where Nk is a normalising factor counting the number of 
wavevectors with the same magnitude, Q m (k), ip m (k) are the 
amplitude and phase of the Fourier component and M m (k) 
are a set of orthonormal basis tensors. Reality of Q is ensured 
by requiring ip m (-k) = -tp m (k) andM m (-k) = M}„(k). 
With these definitions the quadratic part of the free energy 
takes the form 
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The basis tensors M m (k) should be chosen so as to diago- 
nalise F' 2 '. It is convenient to introduce the quantity A k := 



1 — (E.k) and use the direction of the electric field and 
of the wavevector for a given Fourier mode to define a lo- 
cal, right-handed, orthonormal frame field {k, v, w} such that 



E = (E.k)k + A k w. We then make use of the observation 
that the fiexoelectric coupling is antisymmetric to rewrite it 
as +ie fEAkkcadcVd, revealing a formal mathematical simi- 
larity between fiexoelectric and chiral couplings. Combining 
these two coupling terms motivates a transformation to a new 
set of basis vectors, obtained by means of a rotation about w: 
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b := ■±-(e f EA k k + 2q L 1 v) , 
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where /i k = ({2qoLi) 2 + { £ fE At) 2 ) 1 '" is a normalisation 
factor. In this rotated local frame, the choice of basis tensors 
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M = ^(3a<g>a-l) , 
results in the desired diagonalisation of F@\ 
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We see from Eqs. (0]i that there is a strong sense in which 
the transition from purely chiral to purely fiexoelectric cou- 
pling may be viewed geometrically as a tt/2 rotation. Indeed, 
this is already evident in the one dimensional cholesteric and 
splay-bend states of Fig. Q] which are related by a tt/2 rota- 
tion about the vertical direction. This rotation is none other 
than the well known fiexoelectro-optic effect H . 

Given this formal similarity it is natural to ask whether 
a precise correspondence can be made between the stable 
phases, and their properties, in the purely chiral and purely 
fiexoelectric limits. To this end we now focus our attention on 
the purely fiexoelectric sector and analyse a number of struc- 
tures motivated by the analogy with cholesterics. 

The simplest fiexoelectric structure is a one dimensional 
phase that was described in Meyer's original paper [4], which 
we shall call splay-bend, see Fig.QJb). The Q-tensor appro- 
priate to splay-bend is the analogue of that for the cholesteric 
helix in chiral liquid crystals: 
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where we have used the electric field to define the z-axis of 
a Cartesian coordinate system and taken the wavevector to be 
orthogonal to the field so as to maximise A k and hence min- 
imise the free energy. The free energy is readily shown to be 



F = \ (Qi + Ql) - 3Q 2 2 Q h + Q 3 h + (Q2 + Qi) ' 
+ ^(iifc 2 - 2e f Ek) . 
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It is formally identical to that of the cholesteric phase with 
EfE/Li playing the role of the chiral parameter, 2qo- There- 
fore the solutions for Q 2 and as a function of the reduced 
temperature, r, and field strength, E, obtained by minimis- 
ing Eq. ([8]), are the same as those obtained for the cholesteric 
phase [10], provided one substitutes 2qo with SfEjL\. 

Extending the analogy with cholesterics, we now consider 
flexoelectric phases with a higher dimensional structure; flex- 
oelectric analogues of the cholesteric blue phases. Based on 
the heuristic observations that the field picks out a preferred 
direction in space, and that the flexoelectric coupling vanishes 
for wavevectors parallel to the field, it seems likely that the 
free energy will be minimised by two dimensional structures, 
containing only wavevectors orthogonal to the field. There- 
fore we consider first a two dimensional phase possessing 
hexagonal symmetry which may be obtained by choosing the 
fundamental set of Fourier modes to be along the directions 
{±ex, ±(-e a + v / 3e 3/ )/2, ±(-e x - V3e y )/2}. An approxi- 
mate Q-tensor, comprising only this fundamental set, all with 
rn = 2, and a homogeneous component, is given by 
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where H.c. stands for Hermitian conjugate, i/ii,ip2,ip3 € 
[0,27r), 5 <E {0,7r} and the short-hand notation <g> s denotes 
a symmetrised tensor product, i.e., u<g> s := u <g> u. The direc- 
tor field of this hexagonal flexoelectric blue phase is shown 
in Fig. [3] The structure consists of an hexagonal lattice of 
strength —1/2 disclinations separating regions in which the 
distortion is one of pure splay along any straight line pass- 
ing through the centre of the hexagon. Indeed, expanding in 
cylindrical polars (p, <fi, z), the director field in a local neigh- 
bourhood of the centre of each hexagon is 
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and therefore the structure may be aptly referred to as a double 
splay cylinder since it is the clear analogue of double twist 
cylinders in the cholesteric blue phases. 

The free energy of the hexagonal flexoelectric blue phase is 
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This is minimised by choosing the phases to satisfy V'i+'02 + 
ip3=Tr and (5 = 0. As in the case of the splay-bend phase, this 
is identical to the free energy of the cholesteric analogue, the 
planar hexagonal blue phase J3, S ■ Consequently, we may 
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FIG. 3: Director field of the hexagonal flexoelectric blue phase de- 
scribed in the main text obtained from a numerical minimisation at 
paramter values Kb = 1,t — ft! = where the phase is stable, see 
Fig. [4] The structure is periodic in both directions. 



draw on this previous work for cholesterics to conclude that 
the hexagonal flexoelectric phase has lower free energy than 
the one dimensional splay-bend phase at sufficiently large 
field strength. However, in the chiral case the hexagonal blue 
phase is not stable since one of the cubic blue phases always 
has lower free energy. It is therefore necessary to see whether 
a similar result also holds in the flexoelectric case. 

Although one may reasonably expect that the structure of 
the cubic blue phases provides a good starting point for a 
consideration of flexoelectric blue phases, the latter differ in 
a number of respects. In particular, since the flexoelectric 
coupling is vectorial, it defines a preferred direction in space 
which will act to lower the symmetry from cubic to tetrag- 
onal. We therefore considered the set of possible structures 
obtained by selecting the Fourier modes of the Q-tensor to be 

{0, ±%e x ,±jz.e y , ±%(e x + e y ) , ±% (-e x + e,) , 



I 27T 



(-ex + £e,),±£(-e„ + £e,)} , (12) 



where L x , L z are the lattice parameters along the x— and z— 
axes of the conventional unit cell respectively. This set of 
structures allows for considerable scope, containing as spe- 
cial cases Q-tensors that are exact analogues of those repre- 
senting the cholesteric cubic blue phases as well as that of the 
two dimensional square structure shown in Fig.[2jb). We have 
investigated the free energy minima of these structures both 
using approximate analytic calculations and with an exact nu- 
merical minimisation of the free energy for a range of initial 
conditions, varying the ratio of the lattice parameters L x /L z 
and the relative amplitudes and phases of the Fourier compo- 
nents. The numerical minimisation was performed using a lat- 
tice Boltzmann algorithm with an identical technique to that 
described in itldfl . In all instances a single minimum was ob- 
tained which corresponds to the two dimensional square flex- 
oelectric blue phase of Fig.[2jb). 

These results were combined with the analytic expressions 
for the free energy of the splay-bend and hexagonal struc- 
tures, Eqs. ([8]) and (Q~TJ|, to construct an approximate phase 
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FIG. 4: Phase diagram of the fiexoelectric blue phases. The circles 
and solid lines show the numerical phase boundaries, while the ana- 
lytic results are given by the dashed lines. 

diagram. We find that although the two dimensional square 
phase is a local minimum of the free energy within the set of 
structures defined by ( fT2b . it never has lower free energy than 
the hexagonal phase. The same result is obtained numerically, 
leading to the phase diagram shown in Fig. [4] which has been 

plotted in the ke, (t — k|) plane, where ke := 2e 2 ^E 2 / L\, 
in accordance with the conventions adopted in cholesterics 
ifioll . As the electric field strength is increased the hexagonal 
fiexoelectric blue phase becomes stabilised over an increas- 
ing temperature interval between the isotropic and splay-bend 
phases. The discrepancy between the numerical and analytic 
phase boundaries is similar to that found in cholesterics and 
arises from neglecting higher order wavevector harmonics in 
analytic calculations. In the cholesteric blue phases the lat- 
tice periodicity is not identical to the cholesteric pitch, but is 
generally somewhat larger. From the numerical minimisation 
we find the same feature in the hexagonal fiexoelectric blue 
phase, with the lattice parameter being approximately 15% 
larger than the 'pitch' (nLi/efE) of the splay-bend phase at 
the triple point. 

The phase diagram bears a qualitative resemblance to that 
of cholesterics fllOfl . but the details reveal differences between 
chiral and fiexoelectric couplings. In particular, it seems that 
the suppression of wavevectors parallel to the field introduced 
by the quantity Ak is sufficient to convert the fully three di- 
mensional states which are stable at large chiral coupling into 
two dimensional states under fiexoelectric coupling. 

We comment that the local director field in the double splay 
regions of the fiexoelectric blue phases is the same as that in 
the escape configuration lattice phase recently suggested by 
Chakrabarti et. al. However, the mechanism by which the 
structure occurs, and is stabilised, is different. Their proposed 
structure arises without the application of an electric field, be- 
ing instead stabilised by the saddle splay elastic constant and 
including weak anchoring of the director field and surface ten- 
sion at the interface between the nematic and isotropic regions 



of the fluid. Note also that the hexagonal symmetry of the 
phase we obtain is different to the square symmetry predicted 
in that work. 

It is important for potential experimental work to give an es- 
timate of the field strength required to stabilise the two dimen- 
sional hexagonal fiexoelectric blue phase. A good estimate is 
that the strength of the fiexoelectric coupling, etE, should be 
as large as the chiral coupling, 2q^L\, in a compound which 
possesses cholesteric blue phases. This gives a field strength 
E 2qoLi/ef « irK/ep, where we have replaced the 
Landau-de Gennes parameters with their director field equiv- 
alents; an average Frank elastic constant, K, and the average 
of the fiexoelectric coupling constants, e := (e s + eb)/2. p is 
the pitch in the cholesteric phase of a compound which also 
displays blue phases, typically ~ 0.3/xm. Recently materi- 
als have been developed for fiexoelectric purposes and found 
to have a large flexoelastic ratio e/K ~ IV ~ 1 ifxiil . Thus 
in these materials the required field strength is expected to be 
E~ lOVp.m^ 1 , which is within the experimentally accessible 
range. In addition, for the proposed fiexoelectric blue phases 
to be seen, it will be helpful to use a material with as close to 
zero dielectric anisotropy as possible and to take care to en- 
sure that surface effects are negligible compared to the bulk. 
The analogy with cholesterics further suggests that fiexoelec- 
tric blue phases are only likely to be found in a very narrow 
temperature range (~ 1 K) just below the isotropic-nematic 
transition temperature, although this is predicted to expand 
with increasing field strength. 

We thank Davide Marenduzzo for useful discussions. 
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